DESY 97-235 
KAIST-TH 12/97 
BUTP-97/34 
SNUTP 97-169 
December 1997 



Contribution of b ^ sgg through the QCD anomaly 
in exclusive decays B± ^ (r/', 7/)(K±, K*±) and {r]',r]) {K^,K'^] 



A. Ali% J. Chay^ C. Greub"^ and P. Ko'^ 

"Deutsches Elektronen Synchrotron DESY, Hamburg, Germany 
^Department of Physics, Korea University, Seoul 136-701, Korea 
Inst. f. Theor. Physik, Univ. Bern, Bern, Switzerland 
^ Department of Physics, KAIST, Taejon 305-701, Korea 



Abstract 



We compute the decay rates for the exclusive decays {t]' ,'r]){K^, K*^) and 5° — >■ 

{■)]', ri){K^, K*^) in a QCD-improved factorization framework by including the contribution from 
the process b —>■ sgg s{ri',ri) through the QCD anomaly. This method provides an alter- 
native estimate of the contribution b sec s{ri, r]') to these decays as compared to the 
one using the intrinsic charm content of the rj' and 1] mesons determined through the decays 
J/ip ^ [rj,ri' ,7]^)^. The advantage of computing the relevant matrix elements via the QCD 
anomaly governing the transition gg —>■ {r]', rj) is that there is no sign ambiguity in these con- 
tributions relative to the matrix elements from the rest of the operators in the weak effective 
Hamiltonian. Numerically, the QCD anomaly method and the one using the radiative decays 
J/ip ^ [f] , f]' , f]c)'y give similar branching ratios for the decays of interest here. The resulting 
branching ratios are compared with the CLEO data on — >■ rj'K^ and B^ rj'K^ and 
predictions are made for the rest. 



(Submitted to Physics Letters B) 



1. Introduction 



The CLEO collaboration has recently reported measurements in a number of exclusive two- 
body non-leptonic decays of the type B hih2, where hi and /12 are light mesons and the 
inclusive decay 5^ v'^s [fl " 0- particular, large branching ratios into the final states 
including r]' are reported |^: 

B{B^^r]' + Xs) = (6.2± 1.6± 1.3) X 10-^ (for 2.0 GeV <p^' < 2.7 GeV), (1) 
B{B^^7]' + K^) = (6.5lii±0.9) X 10-^ (2) 
5(5° T^'ir") = (4.7+1^ ± 0.9) X 10-^ . (3) 

Interestingly, no decay involving the {riK) or {r],ri')K* modes of either the charged 5^ or the 
neutral 5° (5°) has been observed and the corresponding limits on some of these decays can be 
seen in Of these, the most stringent limit is reported on the decay B^ — > rjK^, for which, 
at 90% C.L., one has § 

B{B^ riK^) < 1.4 X 10"^ (4) 



These measurements have stimulated a lot of theoretical activity - |jT8[ . We will concentrate 
in this paper on the exclusive two-body decays B [t]', rj) {K, K*) for both neutral and charged 
B mesons. 

In ref.. 0, two of us (A. A. and CG.) have studied a number of non-leptonic two-body 
exclusive decay modes of the B^ and B^ mesons in the QCD-improved effective Hamiltonian 
approach, involving the effective four-quark and magnetic moment operators. Since the relevant 
matrix elements of the type {hih2\Oi\B) , where Oi are four-quark operators, are difficult to 



estimate from first principles, one often resorts to the factorization approximation [|T^], in which 
the matrix elements of interest factorize into a product of two relatively more tractable hadronic 
matrix elements. The resulting decay rates depend on a set of effective parameters, which have 
to be determined from experiments. We recall that this generalized factorization approach 
appears to describe the two-body non-leptonic B decays involving the so-called heavy to heavy 
transitions reasonably well [^. Likewise, data on B ^ Kn and B tth decays are well 
accounted for in this framework 0. 

In contrast to the decays B — >■ Kit and B tttt, the decays B {ri,ri'){K, K*) in the 
factorization approach require additionally the knowledge of the matrix elements ('7'| 07^750! 0) 
and (r7|c7^75c|0), emerging from the decay b —>■ s{cc) s{ri',r]). Parameterizing them as 
(^'|c7^75c|0) = —iffj/qfj. and (^|c7^75c|0) = —iffjq^, the quantities of interest for this contribu- 
tion are /^/ and which are often referred to as the charm content of the t]' and r], respectively 



21| . These quantities are a priori unknown but they can be determined in a number of ways, 
also including the 5-decays being discussed here. In [^, these quantities were determined from 
the decays J/ip ^ {rj, rj', rjc)'~i, extending the usual (17, ?7')-mixing formalism p2| to the (?7c, r/', rj) 
system. Using the measured decay widths for the decays J /ip ^ {rj, rj', rjc)'~i and (r^c, rj) 77 
yields \ f^,\ ~ 5.8 MeV and \ f^\ ~ 2.3 MeV 0. 

In the meanwhile, theoretical arguments based on SU (3)-breaking effects in the pseudoscalar 



nonet {n , K , r] , t]') in the chiral perturbation theory approach |2^, and phenomenological anal 



ysis involving in particular the 777 and rj'^y transition form factors |^ , have put to question the 
conventional (one mixing-angle) octet-singlet mixing scheme for the {rj, rj') system. The mod- 
ified two-angle mixing scheme, proposed in p3I], has also implications for B decays involving 



1 



the Tj and rj' meson in the final state. In particular, estimates of |/,^/| and |/^| are expected 
to get revised. Numerically, these quantities depend on the input values of the pseudoscalar 
coupling constants, /o and /§, and the two mixing angles, called 6^0 and 6*8. However, it is found 
that using the best-fit values of the parameters from p4| , which are consistent with Leutwyler's 



estimates of the same ||2^, the estimate of \f^,\ in the modified mixing scheme remains prac- 
tically unaltered. In contrast, the quantity |/^| is considerably reduced due to the small value 
of the mixing angle in the singlet sector, which makes r] an almost octet state. One finds now 
l/^l ~ 0.9 MeV p. With these estimates, it has been argued in that the charm-induced 
contribution h s(cc) s{ri' ,7]) does not dominate the matrix element for rj'K^. The 

resulting branching ratio BR{B^ rj'K^) = (2 — 4) x 10~^ is somewhat lower than but not 
inconsistent with the experimental number in Eq. (2). 

The branching ratio for B^ rj'K^ (and other related decay modes) depends on the sign of 
the quantity f^i (and involving the rj meson), as well as on the phenomenological parameter 
^, which in turn determines the effective Wilson coefficients in the factorization approach 0. 
It is not unreasonable to expect that the value of ^ will be similar in the decays B ^ hih2, as 
the energy released in these decays are comparable; hence this parameter can be determined 
in a number of B decays |2^ . A consistent determination of the parameter ^ will also check 
the consistency of the underlying theoretical framework, namely QCD-improved factorization. 
However, it is desirable to get independent estimates of the quantities f^, and and also settle 
the sign ambiguity present in the method used in |^ . We note that a recent phenomenological 
study has put a bound on /^,, namely —65 MeV < f^, < 15 MeV, with being consistent 
with zero |2^ by analyzing the evolution of the 777 and r]''y form factors, respectively. The 



bounds on from this method are not very stringent and the estimates of this quantity in |^ 
are well within these bounds. 

In this letter we propose another method for computing the contribution of the amplitudes 
b —>■ s{gg) —>■ s{ri',ri). This method is based on calculating the amplitude for the chromomag- 
netic penguin process b —>■ sgg, followed by the transitions gg {1]', rf) which are calculated 
using the QCD anomaly, determining both the sign and magnitude of these contributions. As 
discussed below, the numerical values of /^/ and now depend on the charm quark mass (both 
of them being essentially proportional to Tn'"^). Varying iric in the range 1.3 - 1.5 GeV, we find 
that the QCD-anomaly-method gives f^, = -3.1 (-2.3) MeV and = -1.2 ( -0.9) MeV. 
Hence, in absolute value, turns out to be very close to the one obtained in the {rjc, r]', t])- 
mixing formalism (the two estimates almost coincide for = 1.5 GeV), but the value 
of /^/ is typically a factor 2 smaller in the QCD-anomaly method. The branching ratios for 
B^ {r]',r]){K^,K*^) and 5° ^ {t]' ,ri){K° , K*°) based on the QCD-anomaly method are 
calculated in this letter and compared with the present CLEO measurements and with the ones 
in . We find that the theoretical branching ratios for i?^ —>■ rj'K^ and B^ —>■ rj'K^ are almost 
equal and both are in the range (2 — 4) x 10~^, in agreement with the estimates in [^]. 

2. Estimate of b ^ (r],r]')s via QCD anomaly 

We write the effective Hamiltonian H^q for the AB = 1 hadronic transitions as 
Gp 



/ 10 N 



W=3 



(5) 



2 



where q = d,s. The operators read 
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Here a and /? are the SU{3) color indices and Aj^^, A = 1, 8, are the Gell-Mann matrices. The 
subscripts V ± A represent the chiral projections 1 ± 75. denotes the gluonic field strength 
tensor. The operators Oi and O2 are the current-current operators, 03,...,0e the so-called 
gluonic penguin operators and O7, ...,0io are the electroweak penguin operators. Finally, Og 
is the gluonic dipole operator. In the following we use the next-to-leading logarithmic values 
(NLL) (with respect to QCD) for Ci, ...,Cq, computed in p6|, while the remaining coefficients 
are taken at leading logarithmic precision (LL). Taking as{mz) = 0.118, aewi^z) = 1/128 
and m^op^ = 175 GeV, the coefficients in the naive dimensional renormalization scheme (NDR) 
(evaluated at the renormalization scale n = 2.5 GeV) read: Ci = 1.117, C2 = —0.257, C3 = 
0.017, C4 = -0.044, C5 = 0.011, Ce = -0.056, C7 = -1 x 10~^ Cg = 5 x 10-^ Cg = -0.010, 
Cio = 0.002 and Cg = —0.158. Among the coefficients of the electroweak penguins only Cg has 
a sizable coefficient, which arises mainly from the penguin. 

Working consistently to the precision mentioned above, we include one-loop QCD corrections 
to the partonic matrix elements of the operators Oi, Og and the tree- level diagram associated 
with Og, where the gluon splits into a quark- ant iquark pair. These issues are discussed in detail 
in ref. In addition, we include the photonic penguin diagram associated with the current- 
current operators Oi and O2. All these corrections can be absorbed into effective Wilson 
coefficients C^^ {i = 1, 10) multiplying the four-Fermi operators Oi, Oio in the basis (H). 
While Cf^, Cg^ are given in Eq. (2.5) in ref. 0, the remaining effective coefficients read 



where 



O7 — O7 + -- — Oe , Og — Og , Og — Og + -- — Oe , Oj^q — Oiq , {I ) 



The function AFi is also given in ref. 0. In the following, we often use the following linear 
combinations of effective Wilson coefficients {N^ = 3 is the number of colors): 

a. = Cf + ^Ct^, {z = odd) ; a, = Cf + ^Cf, {i = even) . (9) 
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As described in detail in ref. 0, the hadronic matrix elements {hih2\C^^0i\B) for the two- 
body decays of the form B hih2 are now readily decomposed into various form factors and 
decay constants when applying the factorization approximation. If rj^ ^ is involved in the final 
state, the factorization of and O2 brings in the matrix elements 

(r/^'^|c7M75c|0) (10) 

which have to be estimated. We model them by annihilating the charm- ant icharm pair into two 
gluons, followed by the transition gg ^ rj^^ (see Fig. 1). The first part of this two-step process, 
i.e. b — > s(cc — i> g{ki)g{k2)) which amounts to calculating the charm-quark-loop from which 
two gluons are emitted, has been worked out by Simma and Wyler in the context of a 
calculation in the full theory. Their result is readily translated to our effective theory approach 
and can be compactly written as a new (induced) effective Hamiltonian H^g, 

n I /^eff \ „2 1 

il'I = ~ [Cf + ^) ^V'.K:AM||)^Gf (Z,,G„,). n'd - . (11) 

with G^jy = \ep,yafiG'^^ (^0123 = +!)• In this formula, which holds for on-shell gluons (g^ = 
(fci + ^2)^ = 2fci ■ ^2), the sum over color indices is understood. The function Az5(g^/m^) is 
defined as 

1 r 4 n 2 

Ai5(z) = -1 + - [tt - 2 arctan(- - , for < z < 4 . (12) 

Note that the h — >■ sgg calculation brings in an explicit factor of a^. However, as shown below, 
this explicit factor gets absorbed into the matrix element of the operator resulting from the 
anomaly. So, to the order that we are working, we use the coefficients Cf^ and in eq. ([ll[) . 



By expanding the function Az5(g^/m^) in inverse powers of m^, it is easy to see that the leading 
(1/m^) term in Eq. (11) generates the chromomagnetic analogue of the operator considered by 



Voloshin |2^ to calculate the power (1/m^) correction in the radiative decay B X^^. It 
should be remarked that the corresponding uu contribution in Fig. 1 is suppressed due to the 
unfavourable CKM factors. The tt contribution is included in the effective Hamiltonian via the 
hsgg piece present in the operator Og. However, in the factorization framework, the hsgg term 
in Og does not contribute to the decays discussed. So, the cc contribution in Fig. 1 is the only 
one that survives. 



The (factorizable) contribution from Eq. (|TT]) in the decays B {ri,ri'){K, K*) is of the 



same order (in as) as those of the other operators taken into account in 0. It should be 
remarked that in calculating the amplitude b sgg, there are more contributions in this 
order in «<, than what is shown in Eq. ( [TT| ) and Fig. 1. However, the diagrams where one of 
the gluons is present in the final state yielding b — sg{r],r]') are not expected to contribute 
significantly to the exclusive two-body decays, but rather induce multibody decays. Certainly, 
these configurations have to be included in inclusive decays B — > (r/, ri')Xs [0] but can be 
neglected in the exclusive decays B —>■ {ri,ri'){K, K*). Likewise, configurations in which one 
of the gluons emanates from the effective bsg vertex and the other is bremsstrahled from the 
b or s-quark (or from the spectator anti-quark in the B meson) to form an rj' or rj are non- 
factorizing contributions and they are ignored as the rest of the amplitudes are also calculated 
in the factorization approximation. 

Working out the hadronic matrix element of Eq. (p^ ) using factorization, we now need to 
evaluate the matrix elements 

(r/W|Gf(D^G.,,J|0) . (13) 
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The operator in Eq. ([T3|) can be written as 

We can discard the second term since it is suppressed by an additional power of Qs which follows 
on using the equation of motion, and furthermore, the first term is enhanced by in the large 
Nc limit. The matrix elements of df3{G2^Ga^,a) are related to those of GG; more explicitly 

5/3(^^'^|G'fG„,,jO) = ^(r^(')|G'fG„,,jO) . (15) 

The conversion of the gluons into rj and rj' is described by an amplitude which is fixed by the 
SU{3) symmetry and the axial U{1) current triangle anomaly. The matrix elements for GG 
can be written as 

{v^'^\^GfG^,,M = ml,j;,, . (16) 



In Eqs. (^6]) the decay constants f^, and read 

= ^^^^^^ ' /t/ = sin + cos 6'o , (17) 

where the coupling constants /g, /o and the mixing angles 9s and 9o have been introduced 
earlier. We follow here the two-angle {rj^rj') mixing formalism of ref. [^, where the mass 
eigenstates \ri) and \ri') have the following decompositions: 

\r]) = cos9s\r]s) - sin6'o|?7o), 

\r]') = sin6'8|r78) + cos 6*01770) • (18) 
Collecting the individual steps, the matrix elements in Eqs. ( [T^ ) can be written as 
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(^W(^)|^Gf(D,G„,.)|0)=.g,^/;, . (19) 

One would have naively expected that the gluonic matrix elements are small since they 
contain an extra factor of a^- However, as shown by Eqs. ([19|) , this is obviously not the case 
and the gluon operator with as a whole is responsible for the invariant mass of the rf^^ 
mesons. Also, the combination entering in Eqs. (|19D involving the product of and the gluon 
field operators is independent of the renormalization scale. 



3. Matrix elements for the decays (r^', r7)(K=^, K*=^) and B° {t]',t])(K,K 

To compute the complete amplitude for the exclusive decays, one has to combine the con- 
tribution from the decay b s{cc) — * s{gg) si]^''^ discussed in the previous section with all 
the others arising from the four-quark and chromomagnetic operators, as detailed in 0. The 
resulting amplitudes in the factorization approximation are listed below for all the eight cases 
of interest {rj' ,ri){K^, K*^) and — >■ {rj' ,rj){K^ , K*^). The expressions are given for 

the decays of the B~ and B^ mesons; the ones for the charge conjugate decays are obtained by 
complex conjugating the CKM factors. 
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The coefficients are defined in Eq. @ and the short-hand notation stands for wi 
'U7^(l — 75)^. The quantities and /^/ are given in Eqs. (0), while and /^/ read 

fo 



/,^' = -2^sin^^8 + Acos^o , 



-2 cos^^fi 



sin 9() 



Note that the matrix elements of the pseudoscalar density have been expressed as 



(ii) B- 
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(27) 



It is instructive to compare the matrix elements derived above with the corresponding 
expressions in ||^, obtained by estimating the intrinsic charm quark content in the 77, r]' mesons. 
Concentrating on the decays [t]' ,ri){K^ , K*^), which were the ones worked out in 0, 

and substituting 



c 



(28) 



we get (apart from the small electroweak penguin contributions neglected in @ but included 
above) exactly the same expressions for the decay amplitudes as in [Q. Therefore, we have a 
simple relation between the decay constants f^,, introduced in the intrinsic charm content 
method, and the form factor Ai^ entering via the operator in Eq. (11). The idea of intrinsic 



charm quark content of rj' and r) and the contribution of the operator in Eq. (|ri|) are related since 
this operator comes from the charm quark loop. Using the best-fit values of the [r], 77')-mixing 
parameters from yielding = -22.2°, ^0 = -9.1°, /s = 168 MeV, /o = 157 MeV, which 

in turn yields f^, = 63.6 MeV and = 77.8 MeV, the relations in (|2|) give f^, 3.1 MeV 

(-2.3 MeV) and ~ -1.2 MeV (-0.9 MeV), with m, having the value 1.3 GeV (1.5 GeV). 
The Q CD-anomaly method gives results for which are in good agreement with the ones in 
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(in particular for = 1.5 GeV) but it yields typically a factor 2 smaller value for than 
the method based on the {r]c,r]' ,r])-mixmg Given the uncertainties in the fit parameters 
and approximations in both the methods, our estimates presented here are consistent with the 
parameters f^, and obtained in §]. 

However note that the two approaches are quite different. The advantage of the present 
approach is that the operator in Eq. ([TT|) gives an unambiguous sign relative to the other 
contributions while the method of determining the intrinsic charm quark content via radiative 
decays J/ip ^ {rjc, rj' , 1])'^ can only give the absolute magnitude. In our approach the relative 
signs of the contributions from and O2 to the other contributions are determined; we obtain 
the negative-/^/ (and /^) solution of the two possible ones which were not resolved in p. 

4. Numerical results 

For the numerical analysis, we take the values of the parameters used in Q. The matrix 
elements listed in (i),...,(iv) depend on the effective coefficients Oi, aio, quark masses, various 
form factors, coupling constants and the CKM parameters. In turn, the coefficients and the 
quark masses depend on the renormalization scale fi and the QCD scale parameter A^;^. We 
have fixed A^;^ using as{mz) = 0.118, which is the central value of the present world average 
tts("^z) = 0.118 ± 0.003 The scale fi is varied between fi = nif, and fi = mb/2, but due 
to the inclusion of the NLL expressions, the dependence of the decay rates on /z is small and 
hence not pursued any further. To be definite, we use fi = 2.5 GeV in the following.. The 
CKM matrix will be expressed in terms of the Wolfenstein parameters A, X, p and r]. 
Since the first two are well-determined with A = 0.81 ± 0.06, A = sin 6*0 = 0.2205 ± 0.0018, 
we fix them to their central values. The other two are correlated and are found to lie (at 95% 
C.L.) in the range 0.25 < t] < 0.52 and -0.25 < p < 0.35 from the CKM unitarity fits |32[. 
However, a good part of the negative-p region is now disfavoured by the lower bound on 



the mass mixing ratio AMs/ AMd IQ. Likewise, the ratio -Ri = 0.65 ± 0.40 measured recently 
by the CLEO collaboration §, with Ri = B{B^{B^) n^K^)/l3{B^ n^K^), disfavours 
the region p < Hence, we shall not entertain here the negative-p values and take three 
representative points in the allowed {p,ri) contour with p > 0. These correspond to the three 
values of the CKM matrix element ratio: |Kj6/V"cfc| = 0.08, 0.11, and 0.05, reflecting the present 
central value of this quantity and the upper and lower values resulting from a generous error 
on it. The specific values of p and rj and the legends used in drawing the figures are as follows: 

1. p = 0.05, r] = 0.36, yielding i/p^ + 77^ = 0.36 (drawn as a solid curve) 

2. p = 0.30, r] = 0.42, yielding i/p^ + 77^ = 0.51 (drawn as a dashed curve) 

3. p = 0,?7 = 0.22, yielding p'^ + r]'^ = 0.22 (drawn as a dashed-dotted curve). 

All other curves in Figs. 2-4, through which other parametric dependences are shown, are based 
on using the values p = 0.05, t] = 0.36. The decay constants /^,, /^,, /" and can be obtained 
from /o and /§ by using ^0 and 9^ for the r]'-T] mixing angles. Since q = m\ is rather close to 
the point = 0, and a simple pole model is mostly used to implement the dependence in the 
form factors, we shall neglect the dependence and equate F^^^^q^ = ml) = FQ{*^{q^ = 0). 
The values used for the form factors are listed in Table |I[ 

The quark masses enter our analysis in two different ways. First, they occur in the ampli- 
tudes involving penguin loops. We treat the internal quark masses in these loops as constituent 
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sin 08 1 cos 00 




0.33 


cos 6s sin Bq 
VE V3 




0.28 



Table 1: Form factors at = 0. 

masses rather than current masses. For them we use the following (renormalization scale inde- 
pendent) values: 

mb = 4.88 GeV, = 1.5 GeV, = 0.5 GeV, = m„ = 0.2 GeV. (29) 

Variation in a reasonable range of these parameters does not change the numerical results of 
the branching ratios significantly. The value of rrib above is fixed to be the current quark mass 
value rfTb{n = mb/2) = 4.88 GeV, given below. Second, the quark masses nib, nig, rrid and 
TJiu also appear through the equations of motion when working out the (factorized) hadronic 
matrix elements. In this case, the quark masses should be interpreted as current masses. Using 



mi{mb) = 4.45 GeV |35| and 



m7(l GeV) = 150 MeV , m^(l GeV) = 9.3 MeV , m;:(l GeV) = 5.1 MeV , (30) 



from the corresponding values at the renormalization scale = 2.5 GeV are given in Table 
2, together with other input parameters needed for our analysis. 



rrib 


rus 


rrid 


ruu 


as{mz) 


tb 


4.88 GeV 


122 MeV 


7.6 MeV 


4.2 MeV 


0.118 


1.60 ps 



Table 2: Quark masses and other input parameters. The running masses are given at the 
renormalization scale /i = 2.5 GeV in the MS scheme. 

The branching ratios BR{B^ — >■ rj'K^) and BR{B^ —* rj'K^) are plotted against the 
parameter ^ in Figs. 2 and 3, respectively. Although not indicated by the notation, the branch- 
ing ratios for the neutral i?-meson decays are always understood to be averages with the 
corresponding charged conjugated decays in the following tables and figures. We would like 
to make the following observations concerning the sensitivity of the branching ratios on the 
various input parameters. 

• CKM-parametric dependence: The branching ratio BR{B^ — > rj'K^) shows a mild de- 
pendence (of order 10%) on the CKM parameters whereas BR{B^ —>■ rj'K^) is practically 
independent of them. 

• s-quark mass dependence: The dependence of the branching ratios on the s-quark mass 
is quite marked. To illustrate this we show the two branching ratios calculated using 
m7(2.5 GeV) = 100 MeV through the dotted curves, to be compared with the corre- 
sponding solid curves which are drawn for the same values of the CKM parameters, 
namely p = 0.05 and rj = 0.36, but m7(2.5 GeV) = 122 MeV. The resulting increase in 
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Processes 


BR {i = 0) 


BR (^ = 0.45) 


Experiment 




(2.7 — 3.6) X 10 " 


(2.0 — 2.6) X 10 " 


(6.5::i;4 ± 0-9) X 10 




(2.6 — 5.2) X 10 ' 


(2.0 - 3.9) X 10 ' 


< 1.3 X 10-^ 


B^ rjK^ 


(1.5 — 2.2) X 10 


(0.8 — 1.7) X 10 ° 


< 1.4 X 10 


B± r]K*^ 


(1.5 - 3.0) X 10-6 


(1.6 - 3.1) X 10-6 


< 3.0 X 10-5 


BO _^ ^/j^o 


(3.0-3.7) X 10-^ 


(2.0 - 2.6) X 10-5 


(4.7^^:^±0.9) X 10-5 


BO ^ ^'K*^ 


(2.0 - 6.5) X 10-^ 


(0.8-1.0) X 10-^ 


< 3.9 X 10-5 


BO r]K^ 


(1.2 - 1.7) X 10-6 


(0.8 - 1.1) X 10-6 


< 3.3 X 10-5 


BO ^ rjK*^ 


(2.4-3.4) X 10-6 


(1.9 - 2.8) X 10-6 


< 3.0 X 10-5 



Table 3: Numerical estimates of the branching ratios for the decays B^ —>■ {'r]',r]){K^,K*^) 
and B^ [rj' ,ri){K^,K*^), obtained by varying the CKM parameters in the ranges indicated 
in the text and the s-quark mass in the range 100 MeV < m7(2.5 GeV) < 122 MeV. The 
first column corresponds to using the value of the factorization-model parameter ^ = and the 
second to ^ = 0.45. The third column shows the present measurement and the 90% CL. upper 
limits reported by the CLEO collaboration |[l|, |^ . 

the branching ratios in lowering the value from = 122 MeV to mj = 100 MeV amounts 
to about 20% (and 65% for = 80 MeV). While in literature one comes across even 
smaller values of m7(2.5 GeV), we subscribe to the view that the value m7(2.5 GeV) = 100 
MeV is a realistic lower limit on this quantity. However, if the present high values of the 
branching ratios for B^ — * rj'K^ and B^ — * rj'K'^ continue to persist, one might have to 
consider smaller values of mj. We also expect progress in calculating quark masses on 
the lattice, sharpening the theoretical estimates presented here. 

• Dependence on f^r. Restricting to the negative-/^/ solution, determined by the QCD- 
anomaly method, we plot the branching ratios with = —5.8 MeV, as obtained in 0. 
The other parameters are: m7(2.5 GeV) = 100 MeV, p = 0.05 and ?7=0.35. The results 
are shown through the long-short dashed curves in Figs. 2 and 3. Comparing these curves 
with the corresponding dotted curves, we see that the resulting branching ratios in the 
two approaches are very similar. 

• Dependence on ^: This amounts to between 20% and 35% depending on the other pa- 
rameters if one varies ^ in the range < < 0.5. In all cases, the branching ratios are 
larger for ^ = 0. 

Taking into account the parametric dependences just discussed, we note that the theoretical 
branching ratios BR{B^ rj'K^) and BR{B^ rj'K^) are uncertain by a factor 2. Determin- 
ing the value of ^ from other B ^ hih2 decays in the future, this uncertainty can be reduced 
considerably (see Table H). 

The ratio of the branching ratios BR{B^ — > r]'K^)/BR{B^ — > rj'K^) is a useful quantity, 
as it is practically independent of the form factors and most input parameters. This is shown 
in Fig. 4. We see that the residual uncertainty is due to the CKM-parameter dependence of 



10 



this ratio, which is estimated as about 10%. We get (for < ^ < 0.5) 



0.9-1.02. (31) 



The present experimental value of this ratio as calculated by adding the experimental errors 
in the numerator and denominator in quadrature is 1.38 ± 0.86. Given the large experimental 
error, it is difficult to draw any quantitative conclusions except that the theoretical ratio in 
Eq. (|31|) is in agreement with data. However, we do expect that the experimental value of this 
ratio will asymptote to unity. 

The branching ratio BR{B^ rj'K^) is found to be somewhat lower than the present 
experimental number reported by CLEO. As shown in Fig. 2, we estimate BR{B^ — > r]'K^) = 
(2 — 4) X 10~^ in our framework compared to the experimental measurement of the same (6.5 ± 
1.75) X 10^^. The calculations presented here are in better agreement with the experimental 
measurement of BR{B^ ri'K^), making it to within la. Of course, theoretical rates can 
be augmented by increasing the values of the input form factors given in Table ^ There is 
certainly some room for such enhancement, but in view of the emerging theoretical consensus 
on the estimates of the form factors and the fact that the branching ratios in a number of 
B hih2 decays do not require such enhancement [Q, this can only be modest. Hence, we 
anticipate that the experimental numbers for BR{B^ — > rj'K^) and BR{B^ rj'K^) will 
decrease so as to be more in line with the rest of the CKM-allowed Q CD-penguin-dominated 
two-body B decays and with our estimates! 

We present in Table |] numerical estimates for all the eight branching ratios BR{B^ 
{r]',r]){K^,K*^)) and BR{B^ {r]',r]){K°, K*°)) calculated in the QCD-anomaly method. 
The ranges shown take into account the m^- and CKM-parametric dependence, discussed earlier. 
The entries in column 2 and 3 are based on the choice ^ = 0, corresponding to using Nc = oo 
in the effective coefficients, and ^ = 0.45, which corresponds to the phenomenological value 
estimated in the decays B {D, D*){n, p) [^, respectively 

A number of comments are in order on the entries in Table ^. First, as shown in this table 
and Fig. 4, theoretical estimates of the branching ratios for B^ — ^ r]'K^ and B^ rj'K^ are 
almost equal and they are also the largest for the eight decays shown. So, it is no coincidence 
that these are exactly the ones measured so far. In particular, the branching ratios for the 
decays B^ r]'K*^ and B^ —> rj'K*^ are found to be the smallest in this group, and we 
predict 



BR{B^ r]'K 

BR{B^ 7]'K^) 

BR{B^ 7]'K*^) 

BR{BO T]'KO) 



< 0.02 , (32) 

< 0.03 . (33) 



Likewise, the branching ratios for the decay modes B^ — > rjK^ and B^ rjK^ are smaller 
compared to their //'-counterparts by at least an order of magnitude. We estimate BR{B^ — >■ 
7]K^) = (1 — 2) X 10~^ and a similar value for the neutral B decay mode. On the other 
hand, the branching ratios for the decay modes B^ r]{K^,K*^) and B^ ri{K^,K*^) 
are all comparable to each other and are predicted to be somewhere in the range (1 — 3) x 
10~^. The reason behind this pattern can be seen in the various constructive and destructive 
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interferences involving the eight amphtudes hsted earher. This is in quahtative agreement with 
the anticipation in . 



5. Concluding Remarks 

We have provided an independent estimate of the process b s(cc) s{gg) — > s{ri', rj), us- 
ing QCD anomaly. The resulting branching ratios in B ^ [i]', ri){K, K*) reported here are close 
to the ones obtained in the (r^c, rj' , ?7)-mixing approach [^. The present method also removes the 
intrinsic sign ambiguity inherent in ref . , thereby reducing one source of calculational uncer- 
tainty. Theoretical predictivity is, however, still limited due to the remaining input parameters 
of this framework and we estimate it to be a factor 2. Despite this, very clear patterns emerge 
among the various decay modes considered here, which are drastically different from the ones 
which follow in other scenarios. Hence, ongoing and future experiments will be able to test the 
predictions of the present approach as well as of the competing ones, such as models based on the 
dominance of the intrinsic charm contributions in rj' , as suggested in 0, or models in which 
dominant role is attributed to the soft-gluon-fusion process to form an rj or rj' |T^, |TB[ . In con- 
trast to our approach, these models typically predict similar branching ratios (within a factor 2) 
involving the modes B rj'K and B ^' K* , in both the charged and neutral B decays. Data 
already rules out models predicting BR{B — > rj'K*) > BR{B f]'K), and is tantalizingly close 
to testing also the soft-gluon fusion models which predict BR{B —>■ rj'K*) ~ 0.5BR{B — > rj'K) 
[|T7| , 0. We note that a large intrinsic-charm component in rj' is not substantiated by inde- 
pendent analysis of the 7]'^ form factor ||2^. Soft-gluon-fusion models are not theoretically 



motivated as they show extreme (quartic) sensitivity of the decay widths to the gluon mass 
(an ill-defined quantity) - reflecting that the method being employed in these models is nei- 
ther infrared stable nor predictive. For a definite test of the theoretical framework presented 
here and in 0], more precise data are required and one has to reduce the uncertainty on the 
input parameters, in particular the s-quark mass. However, if in forthcoming experiments, the 
branching ratios presented in Table |^ for the eight decay modes are found to be significantly 
larger (in particular in the modes B^ r]'K*^ and B^ — > rj'K*^), then it will most probably 
be an indication of significant non-factorizing contributions, which may feed into the decays 
B {T]','r]){K, K*) dominantly through the matrix elements of the dipole operator Og. 

In conclusion, we reiterate the two intrinsic assumptions of our approach: (i) absence of 
final state interactions, and (ii) absence of non-perturbative contributions in penguin diagrams, 
which may modify both the magnitudes and phases of the effective coefficients calculated in 
the factorization framework presented here. The first should be a good approximation for the 
decays being considered. For the second, we note that non-perturbative contributions are not 
specific to the decays B — > (r/', ri){K, K*) but are endemic to all such decays where penguins play 
a dominant role Furthermore, we admit that the factorization framework used here and 



elsewhere is vulnerable and it is conceivable that non-perturbative non-factorizing contributions 
do play a significant role in non-leptonic two-body -B-decays being discussed here. This remains 
to be tested experimentally or proven on firm theoretical grounds in a well-defined framework, 
such as lattice QCD. However, it is fair to say that there exists at least a prima facie case in 
some of the related decays, such as i? ^ Kti and B ^ Tm, which support the contention that 
the neglected non-perturbative effects are not overwhelming and the measured decay rates can 
be explained without invoking them 0. Of course, this point of view may have to be revised 
with more precise data. 
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Figure 1: FejTiman diagram contributing to the processes b s{cc) — > s{gg) — > sr]^'^ in the 
full and effective theory. The lower vertex in the diagram on the right is calculated with the 
insertion of the operators and O2 in the effective Hamiltonian approach; the upper vertex 
in both the full and effective theory is determined by the QCD triangle anomaly. 




Figure 2: The branching ratio BR{B^ rj'K^) plotted against the parameter ^. The lower 
three curves correspond to the value m7(2.5 GeV) = 122 MeV and the three choices of the CKM 
parameters: p = 0.05, r] = 0.36 (sohd curve); p = 0.30, r] = 0.42 (dashed curve); p = 0, r/ = 0.22 
(dashed-dotted curve). The upper two curves correspond to the value m7(2.5 GeV) = 100 MeV, 
p = 0.05, ?7 = 0.36 and f^, = —2.3 MeV from the QCD-anomaly method (dotted curve) and 
/^/ = —5.8 MeV from (long-short dashed curve). The horizontal thick solid lines represent 
the present CLEO measurements (with ±lcr errors). 



15 




Figure 3: The branching ratio BR{B° — > rj'K°) plotted against the parameter ^. The leg- 
ends are the same as in Fig. 2. The horizontal thick solid lines represent the present CLEO 
measurements (with ±lcr errors). 
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Figure 4: Ratio of the branching ratios BR[B^ rj'K^)/ BR{B^ rj'K^) plotted against 
the parameter ^. The legends are the same as in Fig. 2. The horizontal thick solid lines 
represent the present CLEO measurements (with ilcr errors estimated as stated in the text). 
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